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ABSTRACT 

Our X-ray study of the nuclear activity in a new sample of six quiescent early-type galaxies, and 
in a larger sample from the literature, confirmed (Soria et al.. Paper I) that the Bondi accretion rate 
of diffuse hot gas is not a good indicator of the supermassive black hole (SMBH) X-ray luminosity; 
in fact, the two quantities appear uncorrelated. Here we suggest that a more reliable estimate of the 
accretion rate must include the gas released by the stellar population inside the sphere of influence of 
the SMBH (generally too small to be probed by Chandra) , in addition to the Bondi inflow of hot gas 
across that surface. We use optical surface-brightness profiles, from archival HST images, to estimate 
the mass-loss rate from stellar winds in the nuclear region: we show that for our sample of galaxies 
it is an order of magnitude higher (typical values of ^ 1O~'*-1O~'^M0 yr~^) than the Bondi inflow 
rate of hot gas, as estimated from X-ray observations (Paper I). Only by taking into account both 
sources of fuel can we constrain the true accretion rate, the accretion efficiency, and the power budget. 
Radiatively efficient accretion is ruled out, for quiescent SMBHs. For typical radiatively inefficient 
flows, the observed X-ray luminosities of the SMBHs imply accretion fractions ~ 1-10% (i.e., ~ 90- 
99% of the available gas does not reach the SMBH) for at least five of our six target galaxies, and 
most of the other galaxies with known SMBH masses. We discuss the conditions for mass conservation 
inside the sphere of influence of a SMBH, so that the total gas injection is balanced by accretion across 
the event horizon plus outflows. We show that a fraction of the total accretion power (mechanical plus 
radiative) would be sufficient to sustain a self-regulating, slow outflow which removes from the nuclear 
region all the gas that does not sink into the BH ( "BH feedback" ) . The rest of the accretion power 
may be carried out in a jet, or advected, depending on the details of the radiative-inefficient accretion 
solutions. We also discuss possible scenarios that would lead to an intermittent nuclear activity, such 
as transitions in the outflow rate or in the efficiency. 

Subject headings: accretion, accretion disks — galaxies: nuclei — galaxies: individual (NGC821, 
NGC3377, NGC4486B, NGC4564, NGC4697, NGC5845) — galaxies: structure 
— X-ray: galaxies 



1. INTRODUCTION 

Most of the galaxies in the nearby universe have in- 
active nuclei (X-ray luminosities < 10**° erg s~^). This 
may be due to an interplay of different factors: low rate 
of gas injection/inflow inside the "sphere of influence"^ of 
the supermassive black hole (SMBH) ; low fraction of the 
available gas being accreted onto the SMBH; low radia- 
tive efficiency of accretion, with the rest of the accretion 
power being advected, or carried out as mechanical lu- 
minosity by a jet or an outflow. 

Our goal is to estimate these factors quantitatively, to 
discriminate between different radiatively-inefficient sce- 
narios, and to outline the power and mass budget inside 
the sphere of influence. To do so, we have selected a sam- 
ple of six quiescent early- type galaxies for which a direct, 
kinematical measurement of their SMBH mass has been 
made (Table 1). Using new Chandra data, we have deter- 
mined the X-ray luminosity of the nuclear sources, the 
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density and temperature of the surrounding hot inter- 
stellar medium (ISM), and the classical Bondi rate (Mb) 
of inflow of the hot ISM into the SMBH sphere of in- 
fluence. The main results are summarized in Table 1, 
see Soria et al. (2005, hereafter Paper I) for details. We 
found that the nuclear sources are much fainter than pre- 
dicted by the standard-disk scenario (which is also ruled 
out by other theoretical considerations at such low ac- 
cretion rates). However, they are brighter than predicted 
by radiatively-inefficient models, perhaps suggesting that 
the Bondi rate is underestimating the true accretion rate 
M (Table 1, Col.9). 

We then considered another eighteen galaxies for which 
the SMBH X-ray luminosity and the hot-ISM Bondi in- 
flow rate have been calculated or constrained from previ- 
ous work (Pellegrini 2005, and references therein; Garcia 
et al. 2005; David et al. 2005). For most of these galax- 
ies, the SMBH X-ray luminosity is lower than predicted 
by the standard ADAF model, suggesting that the true 
accretion rate M ^ Mb- Overall, there is little or no 
correlation between Bondi rate and X-ray luminosity of 
the SMBH (Fig. 14 in Paper I). 

However, the hot. X-ray emitting ISM may represent 
only a small fraction of the gas fuelling the SMBH. That 
can be the case in systems where gas can cool efficiently 
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(cooling timescale < accretion timescale), or, vice versa, 
if gas is injected into the inner regions in a cool or warm 
phase and is accreted before it has time to virialize. Cool 
and warm gas would of course elude an X-ray investiga- 
tion. Furthermore, Chandra observations do not have 
enough spatial resolution to investigate sources of gas 
located inside the sphere of influence of the SMBH (typ- 
ically < 10 pc). In this paper, we use optical brightness 
profiles to obtain a complementary estimate of the gas 
injection rate into the SMBH sphere of influence; from 
the stellar densities and population ages, we estimate 
the characteristic mass injection rates due to stellar out- 
gassing. 

We compare the observed X-ray luminosities of the nu- 
clear sources with those expected from the total accretion 
rates: hence, we constrain the radiative efficiency and the 
fraction of gas that is really accreted, and discuss the pos- 
sible fate of the gas component that does not reach the 
central BH. This allows us to discriminate between differ- 
ent accretion solutions which may be applicable to our 
systems. Finally, we discuss the total accretion power 
budget and the condition for gas mass conservation in- 
side the sphere of influence of the SMBH. We estimate 
the role of feedback-regulated outflows to achieve a mass 
equilibrium inside the sphere of influence, and we show 
that they are consistent with the energy budget of the 
system. 

2. OPTICAL STUDY 

2.1. What we can learn from optical studies 

The main source of hot ISM in elliptical galaxies is 
thcrmalized gas lost by the stellar population, e.g., from 
red-giant winds (Ciotti et al. 1991). Part of the gas that 
may be found in the nuclear region is produced locally; 
an additional contribution may come from a cooling flow, 
from larger radii. Gas can also reach the inner galactic 
regions from the accretion of small satellites (e.g., Pipino 
et al. 2005); this is a fairly common process in ellipticals, 
and is thought to explain the origin of (some) nuclear 
disks. 

Optical studies, with their better spatial resolution, 
can provide a better estimate of the stellar mass, and 
thus the stellar mass losses due to winds, in the circum- 
nuclear region, and therefore a better estimate of the 
gas available to fuel the SMBH. HST data were used 
for this purpose in Fabbiano et al. (2004), to estimate 
the gas available in the nuclear region of NGC821. In- 
ner brightness profiles also allow us to investigate the 
possibility that some nuclei are X-ray faint because of 
a stellar deficit in the core, which may have resulted 
from the gravitational ejection of stars during the coa- 
lesence of two or more SMBHs (Milosavljevic et al. 2002; 
Ravindranath, Ho & Filippenko 2002; Graham 2004). In 
addition, using wider field-of-view ground-based images 
allows us to model the global light-profile and thereby 
acquire an alternative estimate of, and thus check on, 
the SMBH mass using the Mbh^"- relation (Graham et 
al. 2001, 2003; Novak, Faber & Dekel 2005). 

Another reason for studying the optical light-profiles 
is that the optical brightness can be used as a proxy for 
the stellar mass distribution, which in turn, for old pop- 
ulations, is proportional to the low-mass X-ray binary 
(LMXB) distribution (e.g., Kim & Fabbiano 2003, 2004; 
Gilfanov 2004). Therefore, the global optical profiles of- 



fer a useful check on whether the unresolved X-ray emis- 
sion is dominated by faint LMXBs or hot gas at large or 
small radii; and whether the X-ray emitting gas is more 
or less centrally concentrated than the X-ray binaries. 

We have studied the optical light-profiles of the six 
galaxies in our sample. Because of the large apparent 
(optical) extent of these galaxies on the sky, we have 
taken the inner light-profiles from HST images and more 
extended profiles from ground-based images. In the case 
of NGC 4564 we have been able to use the HST hght- 
profilc given in Trujillo et al. (2004) for both purposes. 
Table |21 provides the literature source for each light- 
profile. In general, major-axis i?-band light-profiles have 
been used. 

2.2. Global light-profiles and SMBH masses 

FigureHlshows the best-fitting Sersic (1963, 1968) 
model parameters for each galaxy's light-profile. Sersic's 
/ji/" model is often written as 

/(i?)=/ee-^"[(^/^")''"-i], (1) 
where /c is the (projected) intensity at the (projected) 
effective radius i?c. The term 6„ is not a parameter but a 
function of n chosen in such a way that Re encloses half 
of the total light (Ciotti 1991; Caon et al. 1993; Graham 
& Driver 2005). To minimize the influence of seeing, 
and the spurious effect on the global profiles that might 
have come from possible additional nuclear features or 
depleted stellar cores, we have excluded the inner ~ 1" 
of the profiles from the fitting process, which convolved 
the i?^/" models with the relevant point spread function 
before matching to the observed data. Although gen- 
erally classified as an elliptical galaxy, we concur with 
Trujillo et al.'s (2004, their Appendix C) re-designation 
of NGC 4564 as an SO galaxy. In addition to possessing 
a stellar ring (Trujillo et al. 2004), its light-profile has 
the classic bulge plus outer exponential disk structure 
(Figure 1). 

We mentioned in Paper I that different kinematic 
SMBH mass determinations for NGC 4486B were dis- 
crepant by more than one order of magnitude. More- 
over, Kormendy et al. (1997) report that the use of 
anisotropic models reveals that NGC 4486B may con- 
tain no SMBH at all. To narrow down this uncertainty, 
we decided to use a secondary mass indicator for this 
galaxy. We estimated its mass from the correlation be- 
tween light-profile shape n and nuclear BH mass A/bh- 
Using log(A-fBH) = (3.02 ± 0.47) log(n/3) -f (7.81 ± 0.07), 
with intrinsic variance e = 0.27 dex (Graham et al. 
2005, in prep.), we derive a mass of 0.5tQ;3 x lO^M©, 
as given in Table 2 (Col. 4). We compared this value 
with that expected from the velocity dispersion/SMBH 
mass correlation. For a velocity dispersion a = 170 km 
s~^ (from Hypercat), we expect A/bh ~ 0.6 x 10® Mq if 
we use the correlation in Merritt & Ferrarese (2001), or 
A/bh ~ 0.7 X IO^Mq using the expression in Tremaine 
et al. (2002). Similar values also follow from the bulge 
mass/BH mass correlation. This agreement between all 
indirect indicators suggests that the kinematic mass mea- 
surement obtained from spherical, isotropic models (Ko- 
rmendy et al. 1997) may be an order of magnitude higher 
than the true value. Throughout this paper, we tabulate 
both sets of values (kinematic measurement and Mbh-ti 
value) for NGC4486B. 
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TABLE 1 

Targets of our, study 



Galaxy Type d A/bh log Mb logMEdd log^x log^Edd (M/Mb)a 

(Mpc) (108 M0) (Mq yr-i) (A/q yr-^) (erg s"!) (erg s"!) 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


NGC 821 


E6 


24.1 ± 2.0 


U-O0_Q 35 


-4.42 


0.29 


< 38.7 


46.0 


< 2 


NGC 3377 


E5-6 


11.2 ± 0.5 


1.0«;? 


-4.56 


0.36 


38.5 


46.1 


1.3 


NGC 4486B 


cEO 


16.9 ± 1.3 


[6.0ti:« 


< -3 


1.14 


38.4 


46.9 


> 0.1; 








5+°-^ 
u.a_(, 2 


< -5 


0.06 


38.4 


45.8 


> 2 


NGC 4564 


E6/S0 


15.0 ± 1.2 


n 56+° °=' 


-5.27 


0.11 


38.9 


45.9 


6 


NGC 4697 


E6 


11.7 ±0.8 


1 7+0.2 
-0.1 


-3.82 


0.59 


38.6 


46.3 


0.3 


NGC 5845 


E3 


25.9 ±2.7 




-4.00 


0.74 


39.4 


46.4 


1.4 



Note. — Col.(l): galaxy ID. Col. (2): from the NASA/IPAC Extragalactic Database (NED). Col. (3): from Tonry et al. (2001), except for 
NGC 4486B, from Neilson & Tsvetanov (2000). Col. (4): from Gebhardt et al. (2003), oxeept for NGC 821 (Riehstone et al. 2004), and for NGC 4486B. 
The latter has two alternative values. The first, in square brackets, is from Kormendy et al. (1997), using a two-integral model. A similar result 
(■^■0+4 8^©) obtained with the same method by Kormendy & Gebhardt (2001). The alternative mass determination (a factor of 10 lower) 

was obtained in this paper (Section 2.2), from the Mbh^^ (Sersie index) relation (Graham et al. 2003, 2005 in prep.). Col. (5): Bondi inflow rate 
of the hot ISM into the SMBH sphere of influence; see Paper I for details. Col. (6): Eddington accretion rate, defined in Paper I. Cols. (7), (8): 
unabsorbed X-ray luminosity of the nuclear source, from Chandra observations (Paper I), and Eddington luminosity. Col. (9): accretion rates that 
would produce the observed X-ray luminosities, at ADAF-likc efficiencies, normalised to the Bondi rate. The ADAF model predicts {M / M-b) ~ 0.1. 



In order to check that the secondary mass indicators, in 
particular the MBH~cr and MBH_~n relations, are reliable, 
we have tabulated (Table 2, Cols. 4, 6) the values of the 
SMBH masses estimated for the other five galaxies. In 
all cases, direct and indirect mass determinations agree 
within the errors, which strengthens our confidence in the 
indirect SMBH mass estimates for NGC,4486B, unless 
this is a really peculiar galaxy. 

2.3. Nuclear light-profiles 

We shall now study the stellar distribution in the nu- 
clear region, which will provide an estimate of one source 
of gas available for accretion. The cores of many lumi- 
nous {Mb ^ —20.5 mag) elliptical galaxies are known 
to be partially depleted of stars (e.g., Quillen, Bower & 
Stritzinger 2000; Rest et al. 2001; Ravindranath et al. 
2001; Trujillo et al. 2004, and references therein). In a 
hierarchical universe, giant elliptical galaxies are formed 
from the dry merger of smaller elliptical galaxies (Naab, 
Khochfar & Burkert 2005; Li, Mo & van den Bosch 2005; 
De Lucia et al. 2005). The SMBHs from the pre- merged 
galaxies sink to the center of the newly-wed galaxy where 
they form a binary pair, the orbit of which decays by 
transferring orbital angular momentum to nearby stars 
(Begelman, Blandford, & Rees 1980). Such a gravita- 
tional slingshot evacuates the stars from the center of the 
new galaxy (Milosavljevic & Merritt 2001). A reduced 
number of stars would imply a reduced local source of 
hot gas and may explain the low X-ray luminosities. 

The ffST-resolved, inner light-profiles are shown in 
Figure El Again, the best-fitting Sersic model has been 
shown, although this is only done here to help show devi- 
ations in the data from the relatively smooth outer stellar 
distribution. NGC 4486B is the only galaxy with an ap- 
parent core, rather than a "power-law" profile. We do 
note, however, that even the short (300 s) HST exposure 



of NGC 3377 saturated at the center. A core may there- 
fore exist inside of the innermost usable data point at 
0'.'2. 

The inner arcsec has been excluded from the Sersic 
fit to several profiles. NGC 5845 has a prominent nu- 
clear disk (Section 5.5; Quillen et al. 2000; Trujillo et 
al. 2004) which has been excluded from the fit shown in 
Figure [21 NGC 4697 also contains an additional nuclear 
feature, probably a star cluster, that was excluded from 
the modeling of the underlying host galaxy light, as done 
in Byun et al. (1996). None of the inner light-profiles dis- 
plays evidence of a central deficit of stars. The apparent 
depletion in the core of NGC 4486B is not due to a stellar 
deficit but to a double optical nucleus (Section 2.4). 

Finally, we compared the optical light-profiles and dif- 
fuse X-ray surface brightness profiles (Figure 3) for the 
four galaxies in which we have sufficient counts in the X- 
ray emission. In two cases (NGC 5845 and NGC 4697) 
the optical profiles appear more centrally peaked. A 
detailed comparison and physical interpretation of the 
large-scale (a few kpc) optical and X-ray profiles is be- 
yond the scope of this work. However, we have used the 
X-ray surface brightness profiles in the nuclear region to 
estimate the amount of hot gas available for accretion 
(Paper I); and we shall use the optical light-profiles to 
estimate the rate at which the SMBH can be fuelled by 
stellar winds (Section 2.5). Taken together, these two 
components will give us a good estimate of the total rate 
at which gas reaches the sphere of influence of the SMBH, 
for the six galaxies in our Chandra sample. 

2.4. Peculiarities of NGC4486B 

Among the six galaxies in our sample, one stands out 
for its various peculiarities. NGC 4486B is a rare com- 
pact elliptical (cE) galaxy. The presence of a core in this 
galaxy (Figure 2) contrasts with the cE galaxy M32, 



4 



Soria et al. 



which has a steep inner power-law that continues into 
the resolution limit of the data (Faber et al. 1997, their 
Fig. 2). Although the inner light-profile for NGC4486B 
has come from an F555W (roughly F-band) image, there 
is no sign of dust at its center. The apparent stellar mass 
deficit arises from a lull in flux between two off-centered 
optical nuclei, separated by « 0'.'15 « 12 pc (Lauer et 
al. 1996; Kormendy et al. 1997). Thus, NGC 4486B is the 
only other galaxy known so far with an optical double- 
nucleus around a single SMBH (the other galaxy is M 31; 
Lauer et al. 1993). The two optical nuclei have been in- 
terpreted as the pericenter and apocenter of an eccentric 
stellar disk around the SMBH (Tremaine 1995); such a 
disk might have resulted from the tidal disruption of a 
massive star cluster (Quillen & Hubbard 2003). 

The inner region of NGC 4486B stands out in our sam- 
ple also for its X-ray properties (Paper I). The nuclear 
X-ray source is point-like, with no detection of extended 
structures. However, it is softer than typical AGN (Fig- 
ure 3 and Table 2 of Paper I), inconsistent with a F < 2.5 
power-law spectrum. We do not have enough counts to 
determine whether its X-ray spectrum is a truly steep 
power-law, or a standard AGN spectrum with a soft ex- 
cess (perhaps from dense, hot plasma). As an aside, we 
note that a faint (~ 10^^ erg s~^), soft X-ray source 
was also detected at the optical position of the nucleus 
in the elliptical galaxies NGC 4472 and NGC 4649 (Sol- 
datenkov, Vikhlinin, & Pavlinsky 2003), although it is 
not clear what fraction of that emission comes from a 
central peak in the thermal plasma distribution and how 
much from the accreting SMBHs. And in M31, an (un- 
related?) super-soft source is located at < 1" from the 
nucleus and might easily be confused for the accreting 
SMBH (Garcia et al. 2000) if that galaxy were at the 
distance of NGC4486B. Given its physical peculiarities, 
and the large observational uncertainty in the mass of 
the SMBH, the gas density and the nature of the nuclear 
X-ray source, we arc aware that the accretion proper- 
ties inferred for this galaxy may not be representative of 
typical quiescent elliptical galaxies. 

2.5. Stellar mass-loss rates 

An estimate for the stellar mass-loss contribution to 
the ISM comes from Ciotti et al. (1991) and is given by: 

M.W»1.5xlO-(^) (^^y^M^yr-K 

(2) 

where M*(i) is the rate at which a stellar population of 
age t and B-band luminosity Lb loses gas through stellar 
winds. Mass is also injected into the ISM by Type-la 
SNe, at a rate which depends on specific SN models, and 
is also proportional to Lb, but is only a few percent of 
M, (Ciotti et al. 1991, their Eq. 7). We shah henceforth 
consider this component to be included in M* . In quasi- 
stationary conditions, we expect this gas to flow towards 
the galactic center and to be shock-heated or virialized 
to temperatures kT ~ 0.5 keV, giving rise to the hot 
ISM seen, for example, in five of the six galaxies of our 
sample. For a detailed study of the physical processes 
and timescales involved in the conversion of gas from the 
warm to the hot phase, see Mathews (1990). 
As a first approximation, we shall take the Bondi rate 



Mb (Table 1; definition and discussion in Paper I, Section 
4.3) to represent the rate of inflow of diffuse hot gas into 
the sphere of influence of the SMBH (i.e., across the sur- 
face of radius race); we shall assume that this term takes 
into account the whole stellar contribution from outside 
Tacc- A fraction < 1 of this inflowing gas will then sink 
into the SMBH; the rest may build up inside the nuclear 
region, or be ejected in a wind. In addition, we need to 
estimate how much gas is produced from stellar winds by 
the stellar population directly inside race- 

To obtain an estimate of Lb inside the SMBH sphere 
of influence we need to integrate the optical-luminosity 
density proflles over the sphere of radius race- Firstly, 
we used the luminosity density profiles plotted by Geb- 
hardt et al. (2003) (taking into account the B — V color 
index for each galaxy, from Faber et al. 1997 and the 
Third Reference Catalogue of Bright Galaxies, de Vau- 
couleurs et al. 1991). The only galaxy in our sample 
for which the light-profile is not available in Gebhardt 
et al. (2003) is NGC4486B. We have also performed our 
own deprojection of the light-profiles for all six of our 
target galaxies, using the fitted Scrsic models'"' (see Gra- 
ham & CoUess 1997, their Appendix B). We note that 
a potential complication may arise from our deprojec- 
tion procedure which assumes a spherical stellar distri- 
bution. The presence of flat nuclear disks clearly depart 
from this assumption of sphericity. Nonetheless, our es- 
timated volumetric luminosity Lb inside of race agrees 
with that from Gebhardt et al. (2003) to within a factor 
of two for every galaxy, after correcting for the differ- 
ent color bands. Finally, for each galaxy, we have taken 
the average between the two values of the integrated L b , 
obtained from Gebhardt 's (2003) profiles and our own 
analysis (except for NGC4486B, where only our value 
is available). These values were used in equation (2) to 
determine the mass loss rates. 

Assuming an old population (12 Gyr) for all galaxies, 
we can estimate the mass loss rates within the accretion 
radius of each galaxy (^ 10 pc. Table 7 in Paper I); typ- 
ical values are M* ~ lO'^M© yr"^ (Table 3, Col.3) (see 
also Fabbiano et al. 2004 for NGC 821). In fact, these 
are lower limits: there is evidence that the stellar popu- 
lation in the nuclear region is younger for all six galaxies, 
ranging from « 2.5 to w 10 Gyr (Table 3, Col.2); tak- 
ing the younger age into account yields higher mass loss 
rates (Table 3, Col. 4). As we shall discuss in Section 4.1, 
a fraction < 1 of this gas sinks into the BH, the rest may 
be ejected in a galactic wind, or settle into a nuclear disk, 
or be recycled into other stars. 

In summary, so far we have checked the masses of the 
SMBHs and estimated the rates at which warm and hot 
(Paper I) gas becomes potentially available for accretion, 
inside their spheres of influence. In the next section, we 
shall compare these rates with the rates at which gas 

5 In three of the six galaxies, NGC4486B, NGC 4564, and 
NGC 4697, the Sersic models shown in Figures 1 and 2 provide 
a good fit to the underlying galaxy light, provided that the ad- 
ditional nuclear eomponent in the innermost fa 1" (disk or star 
cluster) is excluded from the fit. While this is the standard pro- 
cedure for the correct determination of n (Table 2), it would lead 
to an under-cstimation of the total flux from inside the sphere of 
influence. So, for the purpose of determining Lg, we refitted new 
Scrsic models to the very inner light-profiles of those three galaxies, 
including also the additional nuclear components. 
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TABLE 2 

Optical surface brightness profiles, Sersic indices and SMBH masses 



Galaxy outer profile n A^bh,!! ""c A^bh ct -'^bh kin inner profile Type Notes 

(10* Mq) (km s-i) (W^Mq) (10* Mq) 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 



NGC 821 


R-band, 1 


4.0 ±0.8 


1 5+2-2 


209 


1 6+°-'^ 
^^■0-0.3 




F160W, 5 


cusp 




NGC 3377 


R-band, 1 


3.0 ±0.6 


Q 7+0.9 


145 


4+°-^ 

"■^-0.1 


l.OtO;? 


F702W, 6 


cusp 


1 


NGC 4486B 


F555W, 2 


2.7 ±0.5 


^+°-^ 

U.O_Q 2 


170 


7+^-2 
-0.2 




F555W, 7 


"dip" 


2 


NGC 4564 


F702W, 3 


2 4+11 


3+1-2 


162 




^6+0-0^ 
'^■^"-0.08 


F702W, 3 


cusp 


3 


NGC 4697 


F475W, 4 


4.0 ±0.8 




177 


o.8l«:« 


-1 7+0.2 

i-'-o.i 


F555W, 7 


cusp 


4,5 


NGC 5845 


R-band, 1 


3.2 ±0.6 




234 


2 r.+0.7 
^■^-0.6 




F702W, 3 


cusp 


1,5 



Note. — Col.(l): galaxy ID. Col.(2): filter and source of the outer light-profile: 1 = Graham et al. (2001); 2 = HST/WFPC2 archive. 
Proposal ID 6099, PI S. Faber; 3 = Trujillo et al. (2004). 4 = HST/ACS archive. Proposal ID 10003, PI C. Sarazin. Col.(3): best-fitting 
Sersic index n for the outer brightness profile (Figure^. Col. (4); predicted SMBH mass from the A/bh — ^ relation (Graham et al. 2003; 
Graham et al. 2005, in prep.). Col. (5): effective velocity dispersion, from Gebhardt et al. (2003), except for NGC4486B, for which we 
take the central velocity dispersion from the Hypercat database. An alternative measurement of the velocity dispersion in NGC 4486B 
is cr fa 185 km s"^ (Kormendy & Gebhardt 2001) which corresponds to a mass of l.Olo 2 ^ 10* Mq. Col.(6): SMBH masses derived 
using the relation in Tremaine et al. (2002). Errors are derived assuming statistical plus systematic uncertainty in the velocity dispersion 
measurements 5% (Gebhardt et al. 2003). Col. (7): SMBH masses derived from kinematic measurements, see Table 1 for references. 
Col. (8): filter and source of the inner brightness profile: 5 = Quillen et al. (2000); 6 = Rest et al. (2001); 7 = Lauer et al. (1995); 3 = 
Trujillo et al. (2004). Col. (9): innermost (resolved) light-profile type. The "dip" in the light-profile of NGC4486B is created by its two 
optical nuclei offset from the central position. Col. (10): 1 = nuclear stellar disk; 2 = double optical nucleus, but a single central SMBH; 3 
= galactic stellar disk; 4 = nuclear star cluster; 5 = dusty disk. 



18 
21 
24 



-0.3 
Au 
0.3 



^ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
r yu.^ n ~ 
?Sw^ 21.49 43.9 4.0 T 

-- NGC M l 

- — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — '■ 


. ' 1 ' 1 ' 1 ' 1 . 

yLtg n ~ 

r in 65.7 4.0 

■- NGC 4697 - 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 _ 

V Me Re ~ 
"X^^ 18.38 2.4 2.7 T 

■- NGC 44865***-««i....,^^^ - 

'■ — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 \ J' T •! 1 — 1 — 1 — - 


r NGC 821 


■ — ^ — \ — ^ — \ — ^ — \ — ^ — \ — ; 

- NGC 4697 - 

> 1 a — H—* — \ — ^ — \ — i 


r NGC 4486B 


r A=0.046 -. 


r A=0.020 -. 


7 A=0.056 ^ *^ ' ^ 


\ — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — ; 

r Me Re n - 
r^'>s..,,,___^ 20.60 44.1 3.0 - 

- NGC 3377 =; 

- — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — '■ 


■ — ^ — \ — ^ — \ — ^ — \ — ^ — \ — : 

V /^o h Me Re n 

18-65 19.5 18.52 7.0 2.4 - 

■ — ^ — \ — ^ — 1 i — t — ^ — \ — ; 


1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — ; 

rV Me Re n " 
7 ^"--....^ 17.80 4.1 3.2 - 

- NGC 5845 

'■ — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — - 


r NGC 3377 


r NGC 4564 


r NGC 5845 

■■i—t—t...i...|.-f"f"f"f"'t"'t-.i..A..*."|- : 


r A=0.034 

- ■ 


7 A=0.039 -. 

1.1,1.1- 


7 A=0.050 -. 

■ .... 1 - 



At 



15 
18 
21 
24 



-0.3 
Au 
0.3 



50 100 
R (arcsec) 



20 



40 60 
R (arcsec) 



80 



10 20 
R (arcsec) 



Fig. 1. — Major-axis light-profiles. The filters used are given in Table|2] The best-fitting Sersic models are shown, while a 
+ exponential disk decomposition is shown for NGC 4564. The residuals of the data about the model are shown beneath each 
used the Sersic indices from these global fits as a check on the SMBH masses (Table 2), using the M^n-n relation (Graham 
2003). 
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actually accretes, indirectly estimated from the observed 
X-ray luminosities. 

3. WHERE DOES THE ACCRETION POWER GO? 
3.1. Accretion rates and efficiencies 

We define a total accretion power Pace = rj'Mc^, where 
T]' is the total accretion efficiency (radiative plus mechan- 
ical), and M the rate at which matter actually accretes 
onto the SMBH. M receives contributions from both the 
hot gas captured at the Bondi radius (Table f and Pa- 
per I), and the additional gas injection from stellar-mass 



losses (Eq.[2] and Table 3). In general, M < (Mb +M») 
because only a fraction of the available gas may reach the 
BH; different accretion models predict different values for 
this fraction. 

Similarly, we define a bolometric luminosity Lboi = 
r/Afc^, where the radiative efficiency rj = fi-rj' < rj' (Paper 
I, Section 5.f). For a standard disk, rj' ^ rj ^ 0.1; for 
purely advective accretion, ry < 77' <C 0.1; in a radiatively- 
inefhcient system dominated by a fast jet or a massive 
outflow, 77' 0.1 and /r <C 1. The X-ray luminosity in 
the 0.3-10 keV band is Lx = /x^boi ~ O.lLboi (Elvis et 
al. 1994; Ho 1999). Finally, we shall use dimensionless 
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Fig. 2. — Major-axis, HST-rcsolved, inner light-profiles. The filters used are given in Table El The long dashed curves show the best- 
fitting (3-parameter) Sersic model. For the lenticular galaxy NGC 4564, an outer exponential disk (short dashed line) has been included. 
For NGC4486B, the best-fitting core-Sersic model (Graham et al. 2003) is also shown (solid curve). Data points excluded from the fit are 
shown by the open circles (see text for discussion). 
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Fig. 3. — Comparison of the optical (curves) and soft X-ray 
(datapoints with error bars) brightness profiles. NGC 4697 and 
NGC 5845 have a relative central spike in the optical brightness 
profile (clearly visible also in Figures 5 and 6 of Paper I) which we 
cannot directly compare with the X-ray thermal-plasma emission, 
due to a lack of spatial resolution. This is one of the reasons why 
we need to take into account both the X-ray and the optical data to 
infer the gas injection rate. In the SO galaxy NGC 4564, the radial 
profile of the diffuse X-ray emission seems to follow more closely 
the optical brightness profile of the virialized stellar component 
(bulge), rather than that of the exponential disk. 



accretion rates to = O.lMc^ / L^dd, "t-b = O-IMbc^ / L^dd 

and TOt = 0.1 (Mb + M^)c^/LEdd- 

We found in Paper I that, for the SMBHs in our six 
target galaxies, Lx/L^dd ^ 10~®-10~'', and that the di- 
mensionless Bondi inflow rates of the hot ISM ttib ~ 10~^ 
(Table 1, and Figure 14 in Paper I). But the hot ISM in- 
flow is only a lower limit to the total gas injection rate. 
From our estimate of the stellar mass losses inside race 
(Section 2.5, and Table 3, Col. 9), we find that the stellar- 
wind contribution is higher than the Bondi inflow rate of 
hot gas for our galaxy sample (all of which have very low 
Bondi rates); this was already pointed out for NGC 821 
in Fabbiano et al. 2004. Hence, we estimate that the di- 
mensionless, total gas injection rate is rh-t ~ 10~^-10~'^ 
(Table 3). 



In Paper I, we compared the observed X-ray luminosi- 
ties of the SMBHs with the inferred Bondi inflow rates, 
and found no clear correlations, as previously noted by 
Pellegrini (2005). In view of the results of our optical 
study, we now suggest that the physical interpretation of 
the faint SMBH emission becomes clearer if we compare 
the X-ray luminosity with the total gas injection instead, 
wherever possible: 

0.1 (Mb + M, 1 ^2 ^ ' vu.iy V"^t 



lO"'"^ - 10""*. 



So, the quiescent SMBH luminosity can be explained by a 
combination of low efficiency and low accretion fraction, 
such that ri{m/mt) ~ 10~^-10~'*. We discuss now which 
accretion solutions may satisfy this constraint. 

3.2. Advective accretion (ADAF models) 

We showed in Paper I that Advection-Dominated Ac- 
cretion Flows (ADAF, Narayan & Yi 1995; Narayan 
2002) under-predict the X-ray luminosities of most 
SMBHs with Bondi infiow rates tob ^ 10"''', but over- 
predict the X-ray luminosities of most SMBHs with 
Tob > 10"''. We shall now compare the observationally- 
determined (unabsorbed) X-ray luminosities with the to- 
tal gas injection rates TOt. Firstly, we recall that in 
the self-similar ADAF model, the luminosity Ladaf ~ 
0.1A/c2(TO/a2) « {m/aYL^dd (Narayan & Yi 1994), 
where the viscosity parameter a ~ 0.1; the radiative 
efficiency 77 ^ O.Ito/q;^ ^ IOto. From a more rigorous 
calculation (Merloni, Heinz & Di Matteo 2003), the X- 
ray luminosity ^x.adaf oc to^- 37^0-97^ jj- ^t^^ ^j^g gg^g jg 
accreted onto the SMBH, even at ADAF-like efficien- 
cies, it would produce X-ray luminosities much higher 
than observed (Table 3, Cols. 6 and 7). Conversely, we 
have estimated the fraction of gas available that has to 
be accreted at ADAF-like efficiencies, to produce the ob- 
served luminosities: we obtain rn/mt ^ 0.01-0.1 for the 
galaxies in our sample (Table 3, Col. 10). 

We have already noted (Paper I) that most SMBHs 
with Bondi rates tob ^ a few 10~ are consistent with 
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radiatively-incfficicnt accretion at a rate m/me ^ 0.1; 
this obviously implies also m/mt < 0.1. We speculate 
that accretion onto this latter group of SMBHs is domi- 
nated by the hot-ISM inflow, i.e., rht ~ me > m* when 
tob ^ a few 10"'*. In the sample of early-type galaxies 
for which luminosity density profiles are available (Geb- 
hardt ct al. 2003, and our own analysis in Section 2.5), 
the integrated B-band luminosities inside the accretion 
radii correspond to mass loss rates m* < a few 10~^, 
for an old or intermediate-age population. Estimating 
the total accretion rate for all the sources in Pellegrini's 
(2005) sample is left to further work. 

We conclude that radiatively-incfficicnt models can ex- 
plain the observed luminosity of our target SMBHs, with 
accretion fraction < 10% of the total gas injected into 
the sphere of influence of the SMBH. Hence, the puzzling 
discrepancy between two groiips of "underluminous" and 
"overluminous" SMBHs (with respect to the ADAF pre- 
dictions) is naturally resolved, when we take into account 
all the sources of fuel. We shall discuss the relation be- 
tween SMBH luminosity and accretion rate in a more 
quantitative way in Section 4, in the framework of a sim- 
ple phenomenological model which will describe the ac- 
cretion power and mass budget. 

In fact, in a few cases the accretion fraction m/mt ^ 
10%. For example, NGC4649 has a dense hot-ISM in 
the central region: rig « 0.5 cm~^, corresponding to a 
Bondi inflow rate tob ^ 10^^'^ (Pellegrini 2005; Sol- 
datenkov et al. 2003). The stellar mass loss rate within 
the sphere of influence is much smaller, m* w 10~^, due 
to the low surface brightness and stellar density in the 
inner region (Gcbhardt ct al. 2003). The observed X-ray 
luminosity of the nuclear source is « 10'^* erg (Sol- 
datenkov et al. 2003); at an ADAF- like efficiency, this 
implies an accretion fraction m/rht ~ 10"'^. For this ob- 
ject, radiatively-inefficient solutions based on convection 
or outfiows (e.g., ADIOS and CDAF models) are strongly 
favoured over the basic ADAF scenario. 

3.3. Advection, convection, outflows or jets? 

Our combined X-ray and optical analysis of quiescent 
SMBHs has suggested that only a fraction < 10% of the 
total gas available from stellar winds and hot-ISM Bondi 
inflow is accreted onto the compact object, and this pro- 
cess is radiatively-inefficient, that is 77 < 0.1. Hence, we 
must now debate what happens to the rest of the energy 
potentially available from accretion, and what happens 
to the gas that does not get accreted. 

In the original ADAF scenario, all the energy that 
is not radiated simply disappears into the SMBH. In 
other varieties of advective flows, for example Convection 
Dominated Accretion Flows (CDAFs: Narayan, Igumen- 
shchev, & Abramowicz 2000), or Advection Dominated 
Inflow-Outflow Solutions (ADIOS: Blandford & Begel- 
man 1999), part of the accretion power is used to sustain 
convection or slow, massive outflows. This scenario has 
the advantage of providing at the same time a physical 
mechanism for reducing the accretion rate, preventing 
most of the gas from reaching the SMBH. W; shall dis- 
cuss the role of outflows in Section 5, when we outline 
the conditions for mass equilibrium inside the sphere of 
influence. 

Another possibility is that most of the non-radiative 
accretion power is carried outwards as mechanical lumi- 



nosity in a fast jet. Models based on a combination of 
radiatively-inefficient accretion flows plus powerful jets 
(Pj ~ O.lMc^) can explain the total energy balance in 
the case of other low-luminosity ellipticals such as M 87 
(Di Matteo et al. 2003), IC4246 (Pellegrini et al. 2003) 
and IC 1459 (Fabbiano et al. 2003). 

It was recently suggested that the dichotomy between 
high-luminosity and low- luminosity AGN corresponds to 
two different accretion modes (Jester 2005, and refer- 
ences therein): thermally dominated (high power from a 
radiatively efficient disk) and non-thermally dominated 
(low power from a radiatively inefficient flow), respec- 
tively. In the latter class of sources, the non-thermal 
emission may come entirely from a steady, compact jet 
(e.g., Falcke, Kording & Markoff 2004 and references 
therein); this state corresponds to the "low/hard" state 
of stellar-mass X-ray binaries. The transition between 
the two modes is expected to take place at « a few 
percent of the Eddington luminosity (Fender, Belloni & 
Gallo 2004), a threshold above which a standard thin 
disk is formed and the jet is quenched. Alternatively, 
this transition could be due to the change in the accretion 
rate caused by an ionization instability (Janiuk, Siemigi- 
nowska, & Szczerba 2004; Siemiginowska et al. 1996). 

Significant support to this scenario has come from 
the "fundamental plane" correlation found by Merloni 
et al. (2003) between X-ray and radio luminosity and 
BH mass in a sample of ~ 100 SMBHs (mostly low- 
luminosity AGN) and a few stellar-mass X-ray binaries. 
None of the six galaxies in our sample have significant 
radio detections; however, five have at least reliable up- 
per limits to their radio-core emission, from VLA obser- 
vations (Section 3.4 in Paper I). Combining the "core" 
X-ray luminosities with the radio flux limits, we note 
that four of the SMBHs are radio faint, falling below 
the best-flt empirical correlation of Merloni et al. (2003) 
(Figure 4). However, the relation has a large scatter 
(ciogfl = 0.88 dex), and our target sample was selected 
to be radio faint. Deeper radio observations are needed 
to determine how discrepant these four SMBHs are from 
those used to derive the fundamental plane correlation, 
and whether or not they belong to a different, radio-quiet 
population without a steady jet, perhaps dominated by 
advection or by slow, massive outflows instead. 

We have also determined the radio-loudness parame- 
ter J?x = i^L^/Lx (Terashima & Wilson 2003), which 
is a ratio of radio core luminosity at 5 GHz over nu- 
clear X-ray luminosity in the 2-10 keV band, useful 
for the classification of AGN and quasars. We obtain 
R < -3.1 (NGC3377), R < -2.6 (NGC4564), R < -3.2 
(NGC4697), and R < -3.3 (NGC5845). These up- 
per limits are consistent with typical ratios expected 
from low-luminosity AGN (Fig. 4 in Terashima & Wil- 
son 2003), extrapolating to luminosities in the ~ 10'^*- 
10^^ erg s~^ range. The conventional boundary between 
radio-loud and radio-quiet nuclei is at i? = —4.5, with 
all low-luminosities AGN being radio-loud. Again, radio 
observations of our targets about one of order of magni- 
tude deeper will determine whether they belong to a new 
kind of X-ray faint, radio quiet SMBH population — or, 
at least, less radio loud than the low- luminosity AGN in 
Terashima & Wilson's (2003) sample. 

The morphologies of the X-ray emission in the nuclear 
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Fig. 4. — Location of five of our target galaxies (those for which 
an upper Hmit to the radio flux Lr is available) in the fundamental 
plane of BH activity The best-fit line and the 90% confidence limit 
are taken from Merloni et al. (2003). The radio flux limits are listed 
in Paper I, Section 3.4. SMBH masses and observed (de-absorbed) 
X-ray luminosities are listed in Table 1. 



regions suggest possible jet-like features in two of our 
target galaxies, NGC 821 and NGC 3377. The former was 
discussed in Fabbiano et al. (2004), who suggest that the 
X-ray flux is consistent with synchrotron emission in a 
jet, or with hot thermal plasma shocked by intermittent 
nuclear outbursts. For NGC 3377, see Paper I, Figure 7. 

4. BALANCING THE MASS BUDGET 
4.1. Disposing of the excess gas 

Since the inflow of gas within the accretion radius 
seems to be much higher than what sinks into the SMBH, 
we want to determine whether the gas builds up until it 
reaches an equilibrium between injection and depletion 
(accretion through the event horizon, plus outflows and 
any other sink terms), and what the equilibrium level 
is; or whether the gas may be depleted in intermittent 
outbursts without reaching a steady state. 

Mass outflows are the most natural mechanism for dis- 
posing of the excess gas and achieving a mass equilib- 
rium (e.g., Di Matteo et al. 2000). We shall determine 
how much mass has to be carried outwards in an outflow 
to ensure a mass balance, how much power is required, 
and what source of energy may power the outflow. Al- 
though we shall focus on outflows for this work, we note 
that scenarios based on star formation in a nuclear disk 
or nuclear star cluster could also in principle provide an 
outlet for the excess gas, without the need for outflows. 
NGC 5845 could be one such example, with its dusty stel- 
lar disk. However, the evolutionary timescale and mass 
deposition rate in nuclear disks are at present unknown 
(Kormendy et al. 2005). 

In the nuclear star-formation scenario, sequestration of 
cold gas into stars could be responsible for the low ac- 
cretion rate onto the SMBH. An observational prediction 
of this model is the presence of a thin disk of dust and 
gas (including cold molecular gas) in the nuclear region, 
with ongoing star-formation; this would be detectable, 
although not necessarily resolved, from the usual indi- 
cators, such as Ha and far-IR emission, and an X-ray 



emission component from SN remnants and high-mass 
X-ray binaries (in addition to the nuclear X-ray source). 
The characteristic size of such disks should be compara- 
ble with the Bondi accretion radius (Paper I, Section 4.3) 
that is < a few 10s pc; their radiative spectrum would 
peak in the mid- or far-IR wavelength range. This sce- 
nario was found to be consistent with the spectral energy 
distribution of M87 (Tan & Blackman 2004), and could 
also be at work in some X-ray faint ellipticals. The far- 
IR, rather than the X-ray, is the most suitable energy 
band for such studies; we have obtained time on Spitzer 
to carry out such an investigation. 

4.2. Outflows and mass equilibrium: a 
phenomenological model 

We shall indicate with Ma{t) the total gas mass inside 
the SMBH sphere of influence (radius race) at a time t; 
with M^{t), the net outflow rate from the same region; 
with Ms, the star plus dust formation rate inside that 
region; with Af, and Mb, the gas injection rate from 
the stellar population and the hot ISM respectively. 
and Mb depend on the stellar and gas content of the 
galaxy: they vary over the Hubble timescale but can 
be taken as approximately constant over timescales < 
10^ yr, i.e., shorter that the evolutionary timescale of 
the stellar population and the cooling time of the X-ray 
emitting gas. For quasi-spherical inflows, the accretion 
timescale < 10^ yr, which is also much shorter than the 
timescale considered here. The total mass Ma is obtained 
by solving the equation 

Ma(t) - Mb + M. - M{t) - AU{t) - AUt), (4) 

where the accretion rate onto the SMBH, M, contains 
contributions from both the inflowing hot ISM and the 
stellar mass losses inside racc^ 

M = aMB + bM,, (5) 

with a < 1 and 6 < 1. In principle, a and b may be 
a function of time and other parameters of the system, 
such as the gas density and temperature. Both from a 
comparison with the observed X-ray luminosities (Ta- 
ble 3) and from theoretical arguments (briefly recalled in 
Paper I), we expect that a 6 < 0.1 for a radiatively- 
inefficient inflow (ADAF, CDAF or ADIOS). Either Mb 
or M* may be the dominant term in the various galax- 
ies, depending on the stellar density and age, and on the 
abundance of hot gas. Henceforth, we shall neglect the 
last term in equation (4) (imposing Ms = 0), i.e., we 
shall assume that the amount of gas deposited in a dusty 
disk or sequestered into new stars (and then new com- 
pact remnants) is negligible compared to the accretion 
and outflow components. 

In the simplest scenario, with constant accretion rate 
(i.e., a and b constant) and constant wind losses, the 
total gas mass inside the nuclear region would increase 
linearly with time unless the outflows are fine-tuned to 
balance the gas sources: 

Mw = (l-a)MB + (l-6)M,. (6) 
To remove the excess gas from the SMBH sphere of in- 
fluence, it has to be imparted a velocity 

v^>v,sc ~ {2GMBBlrf'^ 

«320(r/racc)~^/^(fcr/0.5 keVf'^ km s-\ (7) 
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TABLE 3 

Stellar mass-loss rates, total accretion rates and X-ray luminosities (expected and observed) 



(1) 


A p'p 
(Gyr) 
(2) 


M« 1 
{Mq yr- 
(3) 




(4) 




(5) 


loefO IM* oc^l 
(erg s-^) 
(6) 


InC / /"V A T\ A IT 

lyj^ AJJAb 

(erg s~^) 
(7) 


obs 

(erg s~^) 
(8) 


(9) 


M/Mt 
(10) 


NGC 821 




1.7 X 10" 


-4 


7.1 X 10- 


-4 


3.6 X 10""' 


42.6 


41.3 


< 38.7 


18.7 


< 7% 


NGC 3377 


n 7+0.4 
•^■'-0.5 


5.5 X 10- 


-4 


2.6 X 10- 


-3 


1.1 X 10-3 


43.2 


42.5 


38.5 


92.9 


2% 


NGC 4486B 




[1.3 X 10 


-4 


1.8 X 10- 


-4 


1.2 X 10-s 


42.0 


39.8 


38.4 


(NA) 


5-25%] 






1.4 X 10- 


-5 


1.9 X 10- 


-5 


1.6 X 10-5 


41.0 


39.0 


38.4 


(NA) 


45-55% 


NGC 4564 


^ 8 


1.7 X 10- 


-4 


4.0 X 10- 


-4 


3.1 X 10-'' 


42.4 


41.0 


38.9 


74.1 


12% 


NGC 4697 




5.7 X 10- 


-4 


8.5 X 10- 


-4 


2.2 X 10-* 


42.7 


41.2 


38.6 


5.7 


7% 


NGC 5845 


9 r;+0.4 


4.0 X 10- 


-4 


3.1 X 10- 


-3 


5.7 X 10-4 


43.2 


42.1 


39.4 


31.3 


7% 



Note. — Col.(l): galaxy ID. For NGC4486B, the bracketed line refers to the higher mass estimate (see Table 1). Col. (2): age of the 
stellar population in NGC 821, NGC 3377 and NGC 5845 from Denicolo et al. (2005); age in NGC 4697 from Trager et al. (2000); age in 
NGC 4564 from Sil'chenko (1997); age in NGC4486B from Sanchez-Blazquez (2004). CoI.(3); mass loss rates from Ciotti et al. (1991), 
assuming an age of 12 Gyr. Col.(4): mass loss rates from Ciotti et al. (1991), assuming the ages listed in Col. (2). Col.(5): normalized 
stellar mass loss rate. Col. (6): expected SMBH luminosity if all the gas injected by the stellar population inside the sphere of influence is 
accreted with standard radiative efficiency (r; = 0.1). Col.(7): X-ray luminosity predicted by the ADAF solution, from the model fits of 
Merloni et al. (2003), assuming M = M»_2, and viscosity parameter a = 0.1. Col.(8): observed X-ray luminosity from the nuclear sources, 
from Paper I. Col.(9): comparison between the gas injection rate into the sphere of influence from stellar mass losses and from Bondi inflow 
of the hot ISM. (No diffuse hot gas is detected in NGC4486B, see Paper I.) CoI.(lO): fraction of the (total) injected gas that has to be 
accreted onto the SMBH, to reproduce the observed luminosities at an ADAF-Iike radiative efficiency (i.e., r) ~ lOrh; more precise values 
obtained from the model fit of Merloni et al. 2003). 



spending a power 

For our target galaxies, this is perfectly compatible 
with the energy budget in the system. For example, one 
source of energy readily available comes from Type-la 
SNe, from the same stellar population responsible for the 
gas injection M*. The total SN heating can be estimated 
as 

W)-7.1xl03o.s.(^) (^)"^rgs- 

(9) 

(Ciotti et al. 1991), where ^sn < 1 and s > 1.3. For our 
target galaxies, typical numbers are SN rates ^ 10~^ 
yr"^ inside raco implying heating rates of a few 10'^* erg 
s~^. Prom equations (2), (8) and (9), we obtain that only 
a fraction k « 0.17 x [uw /(500 km s"^)]2 of this power is 
needed to remove « from the SMBH sphere of 
influence. 

Although consistent with the energy available, this sce- 
nario requires an uncomfortable degree of fine-tuning 
between injc;c;tion, outflows and accretion. To make it 
more plausible, we need to make both the accretion and 
outflow rates self-regulated, in order to achieve a quasi- 
steady state. Firstly, to prevent the gas from building 
up indefinitely, we assume that the accretion rate is a 
function of the gas density inside or at the surface of the 



sphere of influence, and therefore, as a first approxima- 
tion, to the total gas mass inside race- We already know 
that this is true for Bondi accretion (M ~ Mb ~ n^, Pa- 
per I; Bondi 1952). If this holds for all gas components, 

we can write M{t) = b'Mg,{t) in Equation (4). The pa- 
rameter {1/b') is the timescalc on which the sphere of 
influence would be drained if all gas sources and outflows 
are turned off. This is of order of the free-fall timescale 
inside the region, ~ rtlcM^^^^ (or a factor of 10 
larger, for ADAF- like inflows: Narayan & Yi 1994), and 
is essentially the same timescale used to estimate the 
Bondi accretion rate. Recalling the physical scaling of 
7*acc (eq. [5] in Paper I), we can write 

where 0.1 < /3 < 1 depending on the details of the ac- 
cretion flow. Secondly, to explain how the wind can be 
fine-tuned to remove ~ 90-99% of the injected gas, we 
assume that the net outflow rate is naturally tied to the 
accretion rate, as described in the next section. 

4.3. Black hole feedback: outflows and jets 

Feedback from the SMBH, either via deposition of me- 
chanical energy or through heating of the ISM, is a pos- 
sible way to stop inflows and/or increase the rate of gas 

outflows from the nuclear region (Ostriker & Ciotti 2005: 
Pellegrini 2005; Omma et al 2004; Ciotti & Ostriker 2001; 
Binney & Tabor 1995). We can model this scenario by 
assuming that a fraction K'i^rf Mc?) = n'rj'b' M.^{f)(? of 
the accretion power is recycled to eject gas from the nu- 
clear region; here r]' is the total efficiency (radiative plus 



10 



Soria et al. 



mechanical) because both kinds of BH feedback may con- 
tribute to the driving of an outflow. Similarly to our 
previous argument, we propose that this feedback power 
will impart a (mass-averaged) velocity > Vesc to the 
outflowing gas. Hence, we require that 



from which we obtain 



(11) 



.500 km s"' ' ^ ^ 

=pb'M,{t), (12) 

where we have taken k' ~ 10% as a characteristic fraction 
of energy feedback, and p is a numerical parameter. By 
inserting A'/„(i) from equation (12) into equation (4), 
and solving, we obtain the asymptotic quasi-stcady-state 
values of mass, average electron density, and accretion 
rate inside the sphere of influence: 



Ma(t = oo); 
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(14) 



M(t = oo) = (Mb + M,)/{p + 1). 



(15) 

In the basic ADAF scenario, p« l/aw 10;p~ 10-10^ in 
CDAF or ADIOS scenarios, which predict a lower mass 
rate reaching the central BH. These values are consistent 
with the observations (Col. 10 in Table 3). 

In terms of the feedback coupling constant k', we ob- 
tain: 

, V. (16) 
500 km s" V 



-) « (^) (-) ( 
0.1 J V7200y \v' J \ 



If this feedback coupling can be achieved (i.e., if k' 1), 
the SMBH feedback suffices to unbind and remove the ex- 
cess gas produced by stellar winds in the nuclear region. 
For a purely advective inflow, {O.lp/r]') « (O.lp/r/) ^ 
a/rh ~ 10^, hence 0.1 ^ k' < 1. In other types 
of radiativcly-inefficient flows, including ADIOS, CDAF 
and jet systems, where a smaller fraction of the accretion 
power is advected, more accretion power is available to 
drive an outflow, so that k' ^ 0.1. 

In our phenomenological model, we require that slow 
(i\v- ^ ^^sc) outflows carry out a fraction p/{p + 1) of 
the available gas mass (not considering a possible mass 
component deposited in a dusty disk or new stars), while 
only a fraction l/(p+ 1) reaches the SMBH. This ratio is 
somewhat reversed if we look at the power budget. The 
same outflows consume only a fraction n' of the avail- 
able total (radiative plus mechanical) accretion power, 
with k' <C 1 for many varieties of radiatively-inefhcient 
solutions. The remaining fraction (1 — k' — /r) of the 
accretion power Pace — r/'Mc^ is still available, for ex- 
ample, for launching a fast jet. 



The kinetic energy carried by a relativistic jet is Pj = 
7jMjc^. a bulk Lorentz factor 3 < 7j < 10 has been ob- 
served in quasars and AGN (Falcke & Biermann 1995; 
Falcke, Malkan & Biermann 1995; Dopita 1997). A 
similar range of Lorentz factors is observed in stellar- 
mass BHs in the "very- high state" (Fender et al. 2004). 
Instead, the range of Lorentz factors for steady jets 
in the "low/hard state" (a case more relevant to our 
sample of quiescent or low-luminosity SMBHs) is still 
hotly debated. Heinz & Merloni (2004) inferred a lower 
limit 7j ^ 5 if the X-ray emission is unbeamed; on 
the other hand. Fender et al. (2004) argued that jets 
in the low/hard state are only mildly relativistic, with 
^'J ^ 0.6c, 7j < 1.4. 

If the power budget of X-ray-faint SMBHs is domi- 
nated by a rclativisitc jet, we expect the mass carried by 
the jet to be « (r7'/7j)M, i.e., only ~ 1% of the mass 
that goes into the BH, and < 0.1% of the mass carried 
away by the slow winds. On the contrary, if the jet is 
slower, for example with uj « 0.5c, its kinetic energy is 
« (l/8)Mjc2, and therefore it can carry a mass Mj ~ M, 
i.e., it can also contribute significantly to the mass out- 
flow. Measuring the radio core luminosity of our target 
galaxies would provide an additional constraint to the 
model. 

4.4. Correlation between X-ray luminosity and accretion 

rate 

In the framework of our phenomenological model, 
we can now also re-interpret the source distribution in 
the X-ray-luminosity versus Bondi-accretion-rate plane 
(cf. Figure 14 in Paper I). We speculate that the appar- 
ent lack of any correlations, noted by Pellegrini (2005), 
and that we confirmed with our additional data, is a 
consequence of the fact that the X-ray luminosity in the 
systems is a function of three physical parameters, rather 
than just one or two, as initially hypothesized. The three 
parameters arc: the hot-ISM inflow rate into the SMBH 
sphere of influence (approximated by the Bondi rate, 
Mb); the stellar mass losses inside the SMBH sphere of 
influence (A'/,); and the fraction l/{p+l) of the total gas 
available that actually accretes onto the SMBH. 

We have replotted the same data points in Figure 5, 
and this time we have considered the dependence on all 
three parameters. We speculate that data points on the 
left-hand side of the plot are characterized by M* » Mb, 
and hence M > Mb and, adopting the self-similar ADAF 
efficiency (Section 3.2), 
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independent of Mb- Data points on the right-hand side 
of the diagram have Mb ^ M*, and hence M/Mb » 
l/{p+l) «; 1 and 
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with < log[Q!(p+ 1)] < 1. We have shown the expected 
location of the nuclear X-ray sources for two choices of 
stellar mass losses (M*/MEdd = 10""^: dashed lines; and 
M^/AfEdd = 10~^: solid lines) and for three possible val- 
ues of the accretion fraction: 1% (indicative of ADIOS or 
CDAF solutions), 10% (ADAF), and 50%. In particular, 
at least five of our six target galaxies have M^/MEdd ~ 
a few 10~^ and are consistent with an accretion fraction 
l/(p-M) ~ 10% (Table 3). 



4.5. Intermittent accretion and state transitions 

So far, we have assumed that the system is in a steady 
state. It is possible to generalize our simple model in 
such a way that the system goes through a cycle of 
low and high luminosity and low and high gas density 
(intermittent accretion, see Section 5.4 in Paper I, and 
Pellegrini 2005), especially if the outflows are powered 
by SMBH feedback. A general condition is that mat- 
ter would alternatively be accumulated inside race (a 
phase in which M + < Mb + M*) and depleted 
(M -h Mw > Mb M*). This can be obtained for ex- 
ample by assuming that the feedback coupling strength 
k' depends on the luminosity or the mass accretion rate, 
being <C 0.01 below a certain accretion rate (thus, unable 
to drive a wind) and > 0.1 above that threshold, with 
the addition of hysteresis or of a time lag with respect to 
the SMBH activity. (For example, the time lag could be 
due to the cooling timescale of the ISM heated/ shocked 
by the SMBH activity.) 

A different way of explaining intermittent episodes of 
nuclear activity is to invoke changes in the relative frac- 
tion of accretion power carried outwards as kinetic en- 
ergy by a relativistic jet (Fender, Gallo, & Jonker 2003). 
We recall that, while most of the mass may be removed 
by a slow wind, most of the accretion power is either 
advected into the SMBH or carried outwards by a fast 
jet. Relativistic electrons in a jet may be responsible 
for inverse- Compton scattering of ambient optical/UV 
photons; the interaction of a jet with the surrounding 
ISM, the accretion inflow, or the slow wind itself will 
produce strong shocks, and associated synchrotron emis- 
sion. Hence, changes in the relative fraction of power in 
the jet will produce noticeable effects in both the radio 
and X-ray emission. 

Another mechanism for intermittent accretion is a 
transition in the physical properties of the accretion 
flow, above a certain density threshold. In fact, it has 
been suggested that changes in the jet power and in 
the inflow structure are physically correlated (Fender et 
al. 2004; Merloni et al. 2003). Siemiginowska, Czerny, & 
Kostyunin (1996) modelled the effect of thermal-viscous 
instabilities in an accretion disk, which can switch be- 
tween a cold state (lower viscosity and M) and a hot, 
fully ionized state (higher viscosity and M). and sug- 
gested that X-ray faint SMBHs may be in the low state. 

Another possible transition is the one between a ra- 
diatively inefficient and a radiatively efffcient solution 
(standard disk, Shakura & Sunyaev 1973). This is ex- 
pected to occur at mass accretion rates m ~ 10~^, at 
which point the ADAF efficiency matches the standard 



efficiency. This corresponds to a mass injection rate 

(Mb + M*) -2.3 X 10-2 (p+l)(^^^) Meyr-i. 

(19) 

If Bondi accretion of the hot ISM is the dominant com- 
ponent, this implies number densities at r = race 

corresponding to hot gas masses > 10 Mq inside the 
sphere of influence of the BH, for uniformely distributed 
gas. 

Quiescent SMBHs may in principle reach this threshold 
at some stage, if the gas inflow rate is as high as in equa- 
tion (19); or even for a lower injection rate, if accretion 
plus outflows plus any other forms of mass sequestration 
cannot balance the gas inflow, so that gas keeps building 
up inside the sphere of influence. At that point, radiative 
cooling inside race becomes efficient, and the gas collapses 
into a standard disk; hence, we may expect the accretion 
rate to increase or jump to a value M ^ (M* — M^) , until 
the disk is drained and the nuclear region is depleted of 
its gas (not only through disk accretion, but also through 
more efficient radiatively-driven outflows). During this 
outburst, the nucleus will appear as a bright AGN, with 
a luminosity Lboi ~ wLEdd ~ 10^'' erg s"^. After the 
disk is emptied, accretion resumes in a low, radiatively- 
incfficient regime, and the gas mass keeps building up 
again until the following outburst. This scenario may 
be applicable to the nuclei of moderately bright systems 
(Lx ~ 10^^-10'''^ erg s"^) but it is almost certainly not 
relevant to our sample of faint elliptical galaxies, and 
probably not for any galaxy of Pellegrini's (2005) sam- 
ple either. Current gas densities and accretion rates in 
our target nuclei arc very far (three orders of magnitude) 
from the transition threshold. 

It is possible that, in some cases, the accretion is out 
of equilibrium and gas may be starting to build up again 
after an outburst has recently emptied the nuclear re- 
gion. If so, applying our phcnomenological model (Sec- 
tion 4.3), we estimate that the mass accretion rate will 
increase, reaching an asymptotic value within < 10® yr. 
However, the observed X-ray luminosities for our target 
nuclei point to accretion rates ~ 1Q~^Mq yr"^, one or 
two orders of magnitude lower than the injection rate. 
It is statistically unlikely that we are observing all our 
targets in the first ~ 10^-10"* yr after each outburst; we 
would, instead, expect to find most of the galaxies close 
to their efficient-accretion threshold or to their asymp- 
totic maximum accretion rate, at any given time. We are 
aware that ours is a biased sample, because we have cho- 
sen some of the X-ray faintest among the nearby galaxies 
with an SMBH mass determination. But we do have es- 
timates of the SMBH luminosity and accretion rate for 
about 1/3 of the « 50 galaxies with a known SMBH 
mass: in almost all of them, the accretion rate is orders 
of magnitude lower than the efhcient-accretion threshold. 
We conclude that they cannot all be out of equilibrium 
or in the initial phase of a transient cycle. 

Even if our quiescent target nuclei could accrete all the 
infiowing hot and warm gas (Table 3) without any out- 
flows, the mass accretion rate would still be at least an 
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Fig. 5. — Normalized X-ray luminosity Lx/i'Edd as a function of normalized Bondi accretion rate me = Me/MEdd (cf- Figure 14 in 
Paper I). Here we have taken into account the presence of stellar mass losses (winds) as an additional source of gas. We have plotted the 
expected luminosities for two "plausible" choices of normalized stellar mass losses inside the accretion radii: Mt/M^dd = 10"* (dashed 
lines) and M,/MEdd = 10~^ (solid lines); this is a typical range of values inferred for galaxies for which the optical brightness profiles 
of the nuclear regions are available (e.g., Table 3; and Fig. 1 in Gebhardt et al. 2003). For each stellar mass loss rate, we have plotted 
the luminosity expected for three different values of 1/(1 +p), which parameterizes the total fraction of gas available that is actually 
accreted onto the SMBH (1%, 10%, or 50%). We have compared these curves with the observed values: Lx and A/b are obtained from 
X-ray observations. Data points in red (in the online version) are the six new galaxies of our study; data points in black are the galaxies 
in Pellegrini's (2005) sample; the data point in green is from David et al. (2005); the large, filled (blue) circles are Sag A* (Baganoff et 
al. 2003) and M31* (Garcia et al. 2005). Assuming a characteristic gas injection rate Mt/M^dd ~ a few 10~*, we conclude that most 
SMBHs are consistent with an accretion fraction ~ 1-10%. However, even lower accretion fractions are inferred for NGC4649; conversely, 
NGC5128 (Cen A) is consistent with an accretion fraction 50%. Another two galaxies (NGC4486B and NGC4261) are consistent with 
accretion fractions ~ 20-50%. 



order of magnitude too low to trigger a transition from 
radiatively inefficient to radiatively efficient flows. Ad- 
ditional sources of gas, such as a more powerful cooling 
flow or satellite accretion from the outer galactic regions, 
or a much younger stellar population in the inner regions, 
would be required for the system to reach the efficient- 
accretion threshhold. For example, a stellar population 
with characteristic age < 500 Myr (equation 2, see also 
Ciotti et al. 1991) would provide a mass-loss rate more 
than an order of magnitude higher, that is probably high 
enough to permit at least transient radiatively-efficient 
accretion. 

5. SUMMARY AND CONCLUSIONS 

We have studied a sample of X-ray faint early-type 
galaxies with Chandra (see Paper I), HST and ground- 
based optical data, with the objective of clarifying the 
relation between X-ray luminosity of their SMBHs, total 
accretion rate, radiative efhciency, and classical Bondi 
rate of capture of the hot ISM by the SMBH sphere of 



influence. We have focused on our new results for a sam- 
ple of six galaxies with kinematic SMBH mass determi- 
nations; in addition, we have used the data available in 
the literature for a larger sample of galaxies (Pellegrini 
et al. 2005; Garcia et al. 2005; David et al. 2005). 

The low SMBH X-ray luminosities rule out radiatively 
efficient, standard accretion. We noted in Paper I that 
radiatively-inefficient accretion models provided better 
estimates, but the X-ray data did not show a clear rela- 
tion between Bondi inflow rates and nuclear luminosities. 
In a few cases, the accretion rates required to match the 
X-ray luminosities are ^ 10% of the Bondi rate; this sup- 
ports radiatively-inefficient flows with convection and/or 
outflows, rather than simple advection. In other cases, 
the required rates are ~ 10% of the Bondi rate, consis- 
tent with basic ADAF models. In a few more cases, and 
in particular for at least five of our six new targets, the 
accretion rates must be > Bondi rate; this is difficult 
to reconcile with the prediction of radiatively inefficient 
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models. 

In this work, wo explained this discrepancy by using a 
new empirical method to estimate the gas injection rate 
into the SMBH sphere of influence. We suggested that 
it can be expressed as the sum of the hot-ISM inflow 
rate (estimated from Chandra observations of the diffuse 
hot gas, with spatial resolutions > 100 pc; see Paper I), 
plus additional contributions from the stellar population 
inside the sphere of influence (stellar winds and Type- 
la SNe). The stellar contribution can be estimated by 
deprojecting the optical brightness proflles to obtain the 
volumetric luminosity densities, and applying standard 
relations between optical luminosity, stellar densities and 
ages, and mass loss rates. We found typical stellar mass 
loss rates --^ 1O~''-1O~'^M0 yr~^; on the other hand, the 
hot gas content varies greatly, leading to X-ray-estimated 
Bondi rates from < 10"^ to ~ IO'^Mq yr ^ over the full 
sample of galaxies. 

We have used these two parameters (Mb and M*) to 
model the total mass injection rate. Only an a priori 
unknown fraction of this gas reaches the SMBH (the 
rest being re-ejected, stored, or turned into new stars), 
which adds another parameter to the model. And only 
an a priori unknown fraction of the accretion power is 
released as X-ray flux (the rest being advected or carried 
out as mechanical luminosity, in a radio jet or a wind). 
Various accretion flow solutions (standard disk, ADAF, 
etc.) have different predictions for the fraction of gas ac- 
creted by the BH, and for its radiative efliciency. Assum- 
ing the ADAF radiative efficiency 77 w lOM/MEdd, the 
observed X-ray luminosities imply that, for most galax- 
ies, only ^ 1-10% of the inflowing gas is accreting onto 
their SMBHs. We suggested that the intrinsic scatter 
and/or observational uncertainty in M^, and in the accre- 
tion fraction is the main reason for the lack of correlation 
between the Bondi rate and the X-ray luminosity of the 
SMBH. 

Measuring the core radio luminosity of an SMBH of- 
fers a possible way of determining what fraction of the 

accretion power is advected and what is instead carried 
outwards as mechanical luminosity by a compact jet, or 
used to power convective flows and outflows (Merloni et 
al. 2003). We only have upper limits to the core radio 
emission of our target galaxies. Most of them lie some- 
what below the Merloni-Heinz-DiMatteo fundamental- 
plane correlation, suggesting that advection, slow out- 
flows, and convection may be more important than rel- 
ativistic jets; however, deeper radio observations would 
be required to provide firm conclusions. 

Based on our X-ray and optical results, we have dis- 
cussed the conditions for mass equilibrium inside the 
SMBH sphere of influence and the fate of the gas that 
does not sink into the SMBH. It is possible that in some 
cases (most notably in NGC5845), part of the excess 
gas may cool down, settle into a dusty/stellar disk, and 
eventually form new stars, even inside the Bondi accre- 
tion radius. In fact, a global understanding of SMBH 
physics (both active and quiescent) requires simultane- 
ous modelling of all three phases (accretion, ejection and 
star formation). However, for this work, we used the 
simplifying assumption that star formation is negligible 
and slow, massive outflows remove ~ 90-99% of the gas 
from the nuclear region. 



Type-la SNe could in principle provide enough power 
to heat the gas and remove it from the central region. 
However, reaching a steady state would require an im- 
plausible degree of fine-tuning between injection, accre- 
tion, and outflow rates. A more likely scenario is that 
both the accretion rate and the outflow rate are self- 
regulating: this happens, for example, if the accretion 
rate is proportional to the gas density or total mass in- 
side the sphere of influence of the SMBH, and the power 
carried by the outflow is proportional to the accretion 
power. In this case, the system can reach an asymptotic 
equilibrium; this mechanism relies on the idea of SMBH 
feedback. We have estimated what fraction of the accre- 
tion power has to be used to drive away the excess gas. 
Pure ADAF solutions require a highly efficient feedback 
coupling, essentially because most of the accretion power 
is lost into the BH. Other radiatively inefficient solutions 
only require that < 1% of the available accretion power 
be used for the feedback. Hence, a fast jet may still carry 
outwards > 99% of the accretion power. If the jet is fully 
relativistic (7J > 5: Heinz & Merloni 2004), it can carry 
outwards ^ 0.1% of the inflowing mass, with the rest be- 
ing either accreted or removed by feedback-driven, slow 
outflows. If the jet is only midly relativistic {v^ < 0.5c: 
Fender et al. 2004), it can carry away an amount of mass 
comparable to what sinks into the SMB, ~ 10% of the 
inflowing gas, with the difference being removed by a 
slower outflow component. As an aside, we note that a 
typical outflow rate of ~ 1O~^M0 yr"'^ from inside the 
sphere of influence would, over a Hubble time, inject a 
gas mass of ~ 1O''M0 into the galaxy: this is negligible 
compared to other sources of gas. 

Alternatively, one can propose scenarios in which the 
system may never reach a steady state: the coupling of 
the central engine with the surrounding ISM may switch 
between an efficient and an inefficient regime, perhaps 
with a time lag or hysteresis, such that gas alternatively 
builds up or gets drained from the SMBH sphere of in- 
fiuence. For example, when active, the SMBH could be 
accreting from the gas accumulated in the inner few pc, 
while at the same time suppressing further inffows from 
larger radii. Observational evidence for feedback-driven 
intermittent accretion in NGC821 was discussed in Fab- 
biano et al. (2004); see also Pellegrini (2005). 

Intermittent accretion may also be triggered by transi- 
tions in the nature of the accreting inffow, regardless of 
SMBH feedback. For example, it could be due to a switch 
from a radiatively inefficient to a radiatively efficient so- 
lution, with the collapse of the gas into an optically thick 
accretion disk; or to a thermal-viscous disk instability. 
Both mechanisms require accretion rates > 10~^MEdd) 
two orders of magnitude higher than what is currently 
observed in our quiescent target galaxies. Instead, they 
may apply to younger elliptical galaxies (< 500 Myr) , or 
to galaxies where most of the accreting gas is provided 
by stronger cooling flows from larger radii, or that have 
undergone recent satellite accretions. 

In this paper, we have only considered the energy and 
mass balance inside the sphere of influence (and neglect- 
ing star formation). A more detailed analysis of the con- 
ditions for SMBH feedback must include a study of the 
timescale on which feedback mechanisms operate (e.g., 
Ostriker & Ciotti 2005; Sazonov et al. 2005). For ex- 
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ample, one has to compare the timescale for Compton 
heating from the central X-ray source with the accretion 
timescale. On the other hand, if the SMBH has a strong 
UV flux component, line-driven winds would be more ef- 
fective to push out the gas originating from stellar mass 
losses (warm phase), before it has time to join the hot 
phase. Shock heating from a jet would also provide ef- 
fective feedback. 

Another related complication we have neglected here 
is that SMBH feedback, in addition to removing gas, 
can also heat up the gas that remains within the Bondi 
sphere, thus making its accretion more difficult (eq.[6] 
in Paper I). In this paper we have simply assumed that 
the diffuse hot gas is approximately isothermal over the 
whole galaxy (see Paper I). If the temperature has a spike 
inside the Bondi radius (not resolved by Chandra) due 
to SMBH feedback, part of the point-like hard X-ray 
emission may in fact be due to the gas itself. In any 
case, it remains true that, in a steady state, the even- 
tual fate of this heated gas is to be either accreted or 
ejected. If it just kept building up, its increased density 
would make cooling more and more efficient (cooling rate 
~ rig). Finally, if nuclear star formation (to be studied 
with Spitzer) is significant, that may also be regulated 
by SMBH feeedback. A study of these issues is left to 
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